Overall Dynamic Properties of 3-D periodic 
elastic composites 

Ankit Srivastava 
Sia Nemat-Nasser 
Department of Mechanical and Aerospace Engineering 
University of California, San Diego 
La Jolla, CA, 92093-0416 USA 

January 26, 2013 

Abstract 

A method for the homogenization of 3-D periodic elastic compos- 
ites is presented. It ahows for the evaluation of the averaged overall 
frequency dependent dynamic material constitutive tensors relating 
the averaged dynamic field variable tensors of velocity, strain, stress, 
and linear momentum. The formulation is based on micromechani- 
cal modeling of a representative unit cell of a composite proposed by 
Nemat-Nasser k Hori (1993), Nemat-Nasser et. al. (1982) and Mura 
(1987) and is the 3-D generalization of the 1-D elastodynamic ho- 
mogenization scheme presented by Nemat-Nasser &: Srivastava (2011) 

We show that for 3-D periodic composites the overall compliance 
(stiffness) tensor is hermitian, irrespective of whether the correspond- 
ing unit cell is geometrically or materially symmetric. Overall mass 
density is shown to be a tensor and, like the overall compliance tensor, 
always hermitian. The average strain and linear momentum tensors 
are, however, coupled and the coupling tensors are shown to be each 
others' hermitian transpose. Finally we present a numerical example 
of a 3-D periodic composite composed of elastic cubes periodically dis- 
tributed in an elastic matrix. The presented results corroborate the 
predictions of the theoretical treatment. 
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1 Introduction 



Recent interest in the character of the overah dynamic properties of 
composites with tailored microstructure necessitates a systematic ho- 
mogenization procedure to express the dynamic response of an elastic 
composite in terms of its average effective compliance and density. 

In the present paper we propose a method of homogenization of 3-D 
microstructured elastic composites which directly provides the overall 
frequency dependent dynamic material parameters. This method is 
the 3-D generalization of the 1-D homogenization scheme presented by 
Nemat-Nasser & Srivastava (2011). The 1-D homogenization scheme 
was shown to be equivalent to the field variable integration scheme of 
Nemat-Nasser et. al. (2010) in the limit of Bloch wave propagation. 
The present method doesn't require the pointwise calculation of the 
values of the field variables within the unit cell. Rather, it gives the 
solution of the elastodynamic field equation to any desired degree of 
accuracy. The method is inspired by the micromechanical homoge- 
nization scheme originally proposed by Nemat-Nasser & Taya (1981) 
and further developed in Nemat-Nasser et. al. (1982). A similar ho- 
mogenization method has been used by Amirkhizi and Nemat-Nasser 
to calculate the effective electromagnetic properties of a periodic com- 
posite Amirkhizi & Nemat-Nasser (2007). 

In what follows we outline our homogenization approach for a gen- 
eral 3-D periodic elastic composite. We show that the resulting av- 
eraged dynamic constitutive parameters are tensorial in nature and 
that the average strain tensor is coupled with the average momen- 
tum tensor. Such a form of the averaged constitutive relation where 
the constitutive parameters (including mass density) are tensors and 
where the average strain (average stress) is coupled with average linear 
momentum has been predicted by Milton &: Willis (2007) and Willis 
(2009). See also Shuvalov ei. al. (2009). Willis has also shown that for 
elastodynamic problems with inelastic strain, the effective properties, 
thus calculated, arc uniquely determined. See Willis (2011). We show 
that this form of the constitutive relation is the natural outcome of 
our homogenization method. Furthermore, we prove, on mathematical 
grounds, that the coupling parameters arising in the averaged consti- 
tutive relations are hermitian transpose of each other, irrespective of 
the material or geometric asymmetries. 
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2 Microstructural homogenization of 
periodic composites 

Here we present a homogenization method based on micromechanical 
consideration of the volume averages of the field variables, viewed as 
measurable macroscopic physical quantities. We express the solution 
to the elasto-dynamic equations of motion as the sum of the volume 
average and a perturbation due to the heterogeneous composition of 
the unit cell, 

Q = Q° + (1) 

where Q represents the field variables, stress (cr) or velocity (ii). The 
aim is to derive a set of constitutive relations for the overall averaged 
parts of the field variables, using the local elasto-dynamic equations 
of motion and constitutive relations. This then provides the homoge- 
nized frequency-dependent material parameters. In what follows, we 
describe our approach using a periodic composite, but the final con- 
stitutive relations also apply to a finite unit cell. 

Consider harmonic waves in an unbounded elastic composite con- 
sisting of a collection of bonded, identical unit cells (17 = {xi : —ai/2 < 
Xi < ai/2; i = 1,2,3}) which repeat themselves in all directions, and 
hence constitute a periodic structure. In view of the periodicity of 
the composite, we have p(x) = p(x + m'l^) and C(x) = C(x + m'l^); 
here m' is an integer, p{x.) is the density, C(x) is the fourth order 
tensor of modulus of elasticity whose inverse is the compliance tensor 
D(x), and I/3,/3 = 1,2,3 denote the three vectors which form a par- 
allelepiped enclosing the periodic unit cell. For time harmonic waves 
with frequency uj, the field quantities are proportional to e^*'^*. For 
waves with wave vector q = gfjXj where Xj is the unit vector in the 
^th (Jirgction and einstein summation convention applies, the Bloch 
representation of the field variables takes the form, 

Q(x, t) = Re [Q(x) exp[i(q ■ x - ojt)]] (2) 

where Q represents the field variables, stress {&), strain (e), momen- 
tum (p) or velocity {ii) whereas Q represents their periodic parts (cr, 
£, p, u). The representation, Eq. [2l separates the time harmonic and 
macroscopic factor from the microscopic part of the field variables. 
Even for a finite unit cell, the Fourier series solution of the microscopic 
part is periodic so that this solution satisfies, Q(x) = Q(x + m'l^). 
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We emphasize here that the frequency and the wavenumber, u and q, 
are, at this point, unrelated and arbitrary. 

The local conservation and kinematic relations are, 

V • cr = -iwp 

(V (g) \x)sym = -ioJE 

where V — >■ V + iq. The corresponding constitutive relations are, 

'" = ^/^ (4) 
p = pu 

where C(x) is the tensor of elasticity and /9(x) is the mass-density. 
These local material parameters represent the structure and composi- 
tion of the unit cell. 

Now we replace the heterogeneous unit cell with a homogeneous 
one having uniform density and stiffness C'^. In order to reproduce 
the strain and momentum of the actual unit cell, field variables eigen- 
stress (5]) and eigenvelocities (U), are introduced. These quantities 
are then calculated, using the basic local field equations and consti- 
tutive relations. The idea stems from the polarization stress or strain 
that was originally proposed by Hashin (1959) and further developed 
by Hashin & Shtriktman (1962,1963) and later by others, in order 
to construct energy-based bounds for the composite's overall elastic 
moduli. The basic tool in these works has been the result, obtained by 
Eshelby (1957) in three dimensions and earlier by Hardiman (1954) in 
two dimensions, that the stress and strain are constant within an el- 
lipsoidal (elliptical in two dimensions) region of an infinitely extended 
uniform elastic medium when that region undergoes a uniform trans- 
formation corresponding to a uniform inelastic strain. 

Here, however, we present a different tool that can be used to ac- 
tually calculate the point- wise values of the elasto-dynamic field vari- 
ables to any desired degree of accuracy. For this, we require that the 
actual values of the field variables at every point within the homoge- 
nized and the original heterogeneous unit cell be exactly the same. To 
ensure this, we require that the following consistency conditions hold 
at every point within the unit cell: 

P = Au-u) 



The eigenstress and eigcnvclocity fields are zero in regions where 
the material properties of the heterogeneous unit cell are equal to the 
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chosen uniform material properties, and p^. From Eqs. ([31 [5]) we 
have, 

V • C° : (V ® <i)sym + oo^P°ii = t^V°U + iuj{\/ ■ S) (6) 

C° : [V O (V • (T)lym + = + iwp°C° : (V ® U),j,™ (7) 

Since the stress and displacement fields (Q) are periodic with the 
unit cell, they can be expanded in a spatial Fourier series, 

Q(x) = qO + = {Q)n + J] Q(^)e^^- (8) 

(Q)n = ^^Q(x)dK (9) 
Q(0 = ^^Q(x)e-^^''dx (10) 

= 8010203 (11) 

^ = ^iXj; = naTr/oo,; n-Q, integers (12) 

where greek indices are not summed. {Q,)n represents the averaged 
value of the field variable over the unit cell and appears in its macro- 
scopic description, and Q(^) represents the local perturbations. 
Eq. (j6l7p become 

- C • C° : (C «) u)sym + uj^P°u = u;V°U - w(C • ^) (13) 

- C" : [C ^ (C • <T)]syn. + OJ^P^CT = - Ujp^'C^ : (C ^ ij)sym (M) 

where C = ^ + For the case of an isotropic reference material we 
have, 

C°H = A°<5ij(5H + p°[5ik6ji + 5u5jk] (15) 

Using the isotropic stiffness tensor, the tensors in Eqs. UflSl [l3|) 
can be inverted to give. 
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u(O = *(O-U(O + 0(O:5](^) 
^(^) = *(^)-U(^)+r(^):S(0 
where (see Appendix), 



(16) 



r2 - r2 



[w^ — c: 



e 



1 



-00 + 



(17) 
(18) 



ijkl 



- 2 



2414-4] 



(19) 



_c2C2][^2 _c2^2- 



■CiCjCkCl + ['5ifc<Jj7 + '^ii'^jA 



VP 



-ujp 



24(4 



-24 



4e 



4e][u^- 



2^21 



4C 



CiCjCp 



(20) 



4e 



where ci = ^yJJ^P~+'2lJP)Jffl is the longitudinal wave speed and ci = 
fjfi / is the shear wave speed. It can be seen that V is symmet- 
ric about its first two and last two indices, $ is symmetric about its 
indices and ^ is symmetric about its first two indices. is not auto- 
matically symmetric about its last two indices but since it contracts 
with eigenstress (Eq. [T6|) which is a symmetric tensor, it can be sym- 
metrized about its last two indices (the antisymmetric part vanishes 
when contracted with the symmetric eigenstress tensor). 



G 



1 



ijp 



(21) 



Now the stress and velocity fields can be expressed as a sum of 
their average and their zero mean periodic components, 



u(x) = {u)+Y,' 



*(0-^ / U(y)e-^ Wy + e(0:i / I](y)e-« Wy 

(22) 
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(23) 

where (u) and (<t) are the average values of the velocity and stress 
fields, respectively, taken over a unit cell. 

To make the homogenized unit cell point-wise equivalent to the 
original heterogeneous unit cell, the homogenizing fields are required 
to satisfy the following consistency conditions: 

D(x) : [{a) + aP] = 0° : [{a) + o-^' - S] (24) 



a(x) = (^)+^e^«- 



p(x)[(u)+uP]=/,0[(u)+u^'-U] (25) 

where D and D'^ are the compliance tensors of the actual and the 
reference materials respectively. The periodic parts of the velocity 
and stress fields, from Eqs. (f23| [22]) . are now substituted into the 
above equations. This gives a set of 2 coupled integral equations 
which yields the required homogenizing stress and velocity fields that 
exactly and fully replace the heterogeneity in the original medium. 
Our immediate concern, however, is not the point-wise representation 
of the heterogeneous medium, but, rather, it is the determination of 
the averaged field values, though the solution technique also yields the 
point-wise values of the field values as well. 

We now divide the unit cell into a number of subregions, 0^. Then, 
we average the perturbation fields over each such subregion to obtain. 



{a^n^ = ^P» = ^ f aP{^)dV. (26) 



(u^)f.. =vP^ = ^ I vP{^)dV^ (27) 
= Y^g^ii) • ^ / U(y)e-^«-Wy + e(0 : 1 / S(y)e-^ydyy 
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= 7^ / e^^ 'dy, (28) 

We now replace the integrals in Eqs. (pHj [27|) by their equivalent 
finite sums and set, 



Eqs. (1261 [27l) then yield the following expressions: 



(29) 



(30) 



(31) 



where the repeated index, /3, is summed over the number of subregions, 
/3 = 1, . . . , a, and greek indices serve as labels for tensors rather than 
components of a particular tensor. The coefficient tensors in the above 
equations are defined by, 

M"/3 = ^5"(^)/V(-OM(0 (32) 

In these equations, (3 is not summed. Averaging the consistency 
conditions over each subregion a and using Eqs. (130} [3T]l . we have. 



r(u) 

~a/3 



/"^M"^ ; a not summed 

p° — 



(33) 



For each tensor in Eq. (j33p we have, 
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M"^ = Yl [r5"(^)/V(-^)M(^) + /"g"(-^)/V(^)M(-^) 

= E [/V(0/V(-OM(0 + /V(-0/V(OM(-0 



(34) 

~ a/3 

where * denotes conjugation. This property also holds for tensors T 
and since they are derived from f"^ and by modifying only 
those tensors for which a = /3. We also have the following identity 
(see Appendix), 

^ mnp ^mnqr^qrp r ^pmn ^pmn \'~'^J 

The above identity suggests that it is expedient to express Eq. (j33|) 



in terms of the modified tensors, ^ and 0. This can be achieved by 
multiplying the matrix corresponding to the tensor vE' with the stiffness 
matrix and multiplying the matrix corresponding to the tensor 
with a diagonal matrix consisting of on its diagonals. We also 
transform tensor T in a similar way in order to yield T which in 
addition to having minor symmetries over its first two and last two 
indices, also has major symmetry over them (See Appendix), 

^mnkl = ^mnij^ijkl (36) 

and transforming $ to 

= P^hj (37) 



To facilitate inversion and solution, Eqs. ()33p are now expressed 
in their equivalent matrix form, 
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[fi]{(^)} = -[w]-^ [n 
[f^]{(u)} = -[/>' 



[W] {S} - [W]"' [D° 



{U} 



3(q- 
3(a- 



1)+1 
l)+3 

1)+1 







{U} 



rv\6("-l)+6 
l-^-l6(a-l)+l 
3(a-l)+3 



f2 



= {(^22) (0-33) (0-23) (cr3l) (0-12)}'^ 

- {{ni) {U2) (na)}^ 

fy-ia y-iQ; yia v^a y^Q v-iCk \T 
1^11 ^22 ^33 ^23 ^31 ^12/ 

6(q - 1) + 6 
6(q - 1) + 6 
3(a - 1) + 3 
3(a - 1) + 3 



{Uj'3(a-l)+l 
f.l. — f"A. •■ 

i = 6(a - 1) + 1 
j = 6(q - 1) + 1 
= 3(a - 1) + 1 
/ = 3(a - 1) + 1 
Q = 1 : a 

(38) 

where [D''] is a 6a x 6a symmetric matrix with the 6x6 comphance 
matrix [D''] about its diagonal and [W] is a 6a x 6a symmetric ma- 
trix with the 6x6 diagonal matrix [W]=diag(l 1 1 2 2 2) about 
its diagonal. Similarly, [p'^] = [l/p^] is a 3a x 3a symmetric ma- 
trix with the reciprocal of the reference material's density, , on its 
diagonal. 



In the above equations 



is a 6a X 6a square matrix and 



a 3a X 3a square matrix. For any a and /3 we have, 



IS 



■pop 
pqrs 

- - 
PI 



1^ rspqS 



X'- pqrs J 



(39) 



where * denotes conjugation. The above result is general since for 
a = /3 tensors T and $ are purely real and symmetric so that, 
combined with Eq. ([39]), the corresponding square matrices are her- 
mitian. 

Using Eq. (|35p we also note the following identity for the matrices 



and 







^ mnp 



9^/3 
^pmn 







(40) 



More generally, denoting a hermitian transpose by f, we have the 
following identities for the matrices. 
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r 





(41) 





t 












To expedite notation we omit the square brackets denoting a matix 
in further calculations. Curly braces and angle brackets indicating 
vectors and averaged vectors resectively are retained. We express the 
solution to Eq. (j38p in the following form: 



{S} = A(^) + A(u) 
{ij} = S{a) + n{u) 



where, 



(42) 



A 

A : 



W-^(D°)"^f W + W"^(D°)"^ 
-W"^(D°)"^fW+ 



0W 



W i(D") 0W 



"^W-1(D0)-i** ^pOf2 (43) 



* + (p°)^0r * {p")-'&T n"wi 



-1 



-1, 



* + (p°)^0r * 



-1 



Eqs. (1421) expresses the vectors of eigenstress and eigenvelocity in 
each subdivision in terms of the vector of average stress and average 
velocity. It can be seen from Eq. (I33p that the vectors of the average 
quantities in Eqs. ([I2|) are composed of a times repeated vectors 
of the averaged quantities. Therefore, Eqs. (|42]l can be condensed to 
express the vectors of eigenstress and eigenvelocity in each subdivision 
in terms a 6 x 1 vector of average stress and a 3 x 1 vector of average 
velocity. Furthermore, we can also average the vectors of eigenstress 
and eigenvelocity over all the subregions to finally get. 



{<S)}6xl = (A)6x6(cr)6xl + (A)6x3(u)3xl 
{(U)}3xl = (S)3x6(cr)6xl + (J^)3x3 Wsx 1 



(44) 



where. 
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(A) 
(A) 



-1 ^ 



1l^T 



W"^(D^V rW + W"^(D")-^** 0W 



2\T 



(H) = -(F^) 



^ - -1 c 



(0) = (F2f -(pO)-i* + (pO)-i0r ^* 

where, 



-1 



(45) 



[F^] 



[F2 



/I 

0/^ 
0/^000 

f 



/loo 

0/^0 
0/1 
f 



Now we average the consistency conditions over the entire unit cell 
to express the average strain and average momentum in terms of the 
average stress and average velocity tensors. Noting that the average 
of the periodic parts vanish when taken over the entire unit cell, we 
have, 



(£)=DOw[(a)-(I])] = D(cr) + Si(u) 

(p)=p°[(u)-(U)] = S2(cr) + p(u) 
The effective parameters are given by, 

p = D°W [1 - (A)] 
Si = -D°W(A) 
S2W = -pO(H) 

p = po[i - m 



(46) 



(47) 



(48) 
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Eqs. (j461 H7|) are our final constitutive relations for the homoge- 
nized composite. These equations can be expressed in the following 
tensorial form, 



(e) = D : {a) + Si • (li) 



(49) 



(p)=S2:M + p-(u) 



(50) 



It is shown in the following section that the effective compliance 
tensor D is such that D^jf^i = Dl^^ - and the effective density tensor p 



the relation {82)*^]^ = {Si)jki- The effective constitutive relation can 
alternatively be expressed as (See Willis 2009, 2011), 



where C = D"!, S = -C : Si, S = S2 : C, and = p - S2 : C : Si. 
Willis has shown that the constitutive relation expressed in the above 
form has a self-adjoint structure, i.e., S(q,a;) = S(— q, (x)). Using the 
properties of the effective parameters proved in the next section, it can 
be shown that the self-adjointness of the structure of the constitutive 
relations of Eqs. (j50|52p emerges identically from our formulation. 

2.1 Properties of the Constitutive Parameters 

We note that the 6x6 matrix D'^ in Eq. (j48p can be taken inside the 
matrix multiplication in Eq. (j45|) by converting it to a 60; x 6a matrix 

which is composed of D'^ about its diagonal. Similarly W is taken 
inside the matrices by converting it to W and is taken inside the 
matrices by converting it to . 

We have, 




(^) = C : (e) + S • (u) 



(51) 



(p) = S:(s)+pi.(u) 



(52) 
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D = D°W [1 - (A) 



D°W - D°W(A) 



D°W - D^WfF 



1a 



0W 



il\T 



w-i(D°)-ir(D' 



Q\-l 



1 A 



GfD 



0\-l 



-1 A 



D'JW-(F^)^ -(D")-^r(D^O" +(D")"'** 0(D^) 



-1 



Considering Eqs. (|41|) and the fact that D" and W are hermitian, 
the above treatment shows that D is also hermitian. It can be shown 
that for a hermitian matrix M, transformations of the form (F^)"'"MF^ 
or (F^)^MF^ result in hermitian matrices. This combined with the 
fact that D*^W is real and symmetric proves that D is also hermitian. 
Similarly it can be shown that p is also a hermitian matrix. 

Taking a hermitian transpose of Eq. (148 and multiplying both 
sides by W"-*^ on the left we have, 



-W-i(H)tpO 



W 



-1 



(f2)V -(p")"'* + (p'')"'er (p'J)-^0r^D"F 



(F^f r-(pO)-i* + (pO)-i0f \p^)-^&t ^d°wfi 

1-1-1 



^Ot:.1 



re * + * 



+ 



-1 



P°f2 



t 



and, 
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-D°W(A) 



+W-i(D°)-i\E'# 



^0W 



-1 



rw + 0w 



0t2 



r + 



-1 



(54) 



proving that = Si. Since the matrices [D] and [p] are hermitian, 
we have Dij^i = D'^nj and pij = p*^ where * denotes a complex con- 
jugate. This means that the scalar given by a = {(^)ijD*jj^i{cr)'^i = 
{a)ijDkiij{a)li has a complex conjugate a* = {'y)kiDlnj{a)*j. Since 
the pairs i,j and k,l are symbols upon which summation is carried 
out, we find that a = a* or that a is real. Similarly {u)ip*j(u)*j can 

be shown to be a real scalar. The relation = Si implies that 



ijk 



{Si)jki- Now energy which can be written as, 



E 



(55) 



2 [(,'^)i3^*jkl{'^)kl + '\'^)ij{Sl)*jkiu)l + {u)i{S2)*jk{'^)*k + {u)ip*j{u)*\ 

\ [{(^)ijD*jki{(^)*ki + {cr)ij{Si)*jk{u)*k + {u)i{Si)jki{cr)*k + {u)ip*j{u)*] 
I [{a)i,D*^,i{a)li + [{a),,{S,)*^,{u)l + ((cT),,(5i)^,(n)^)*] + {u)iP*^{u)* 



(56) 



is proved to be real-valued. Taking into account the symmetries of the 
homogenized constitutive parameters, it must be noted that the aver- 
aged constitutive relations of Eqs. (j46l I47p can have 45 independent 
constants at the maximum. Depending upon the material properties 
of the constituents and the symmetry of the unit cell, the number of 
independent constants will vary. The essential relations of the consti- 
tutive tensors proved in this section will hold regardless of the material 
anisotropy or directional asymmetry. 
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3 Numerical Example: Elastic cubes 
periodically distributed in an elastic ma- 
trix 

We consider a composite with a 3-D periodic microstructure composed 
of the periodically distributed unit cell shown in Fig. ([T|). 



1 cm cube 




□ Young's Modulus: 2x10' Gpa □ Young's Modulus: 200x10' Gpa 
Density: 1000 kg/m^ Density: 3000 kg/m^ 

Poisson's ratio: .3 Poisson's ratio: .3 



Figure 1: Schematic of a 3-D periodic composite 

The composite is periodic in all three directions and the wave is 
assumed to propagate in the positive x direction. Due to the peri- 
odicity of the composite, Bloch waves with wavenumber q in the x 
direction propagate in the structure. The wavenumber is related to 
the frequency by the dispersion relation. 
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Normalized Wavenumber (Q) 

ure 2: Dispersion curve showing the first five modes of the 3-D composite 

Fig. ([2]) shows the dispersion results for the first 5 modes calculated 
from the mixed variational formulation (See Iwakuma & Nemat-Nasser 
1982, Nemat-Nasser &: Srivastava 2011, and Nemat-Nasser et. al. 
2011). and indicate positive shear modes and B^ indicates 
positive longitudinal mode. B^ and B^ are shear modes for which 
the group velocity is anti-parallel to the phase velocity and are called 
negative shear modes. B^ and B^ (B^ and B^) lie on top of each other 
because due to the symmetry of the unit cell, the shear directions, y 
and z, are indistinguishable for a wave propagating in the x direction. 

Table ([3]) shows the effective properties calculated from the mi- 
cromechanical homogenization method proposed above. The effective 
properties are calculated for Q = 1.5 on the positive longitudinal 
branch of the dispersion curve (^3). The results were calculated by 
taking the reference material to be the matrix so that discretization of 
the inclusion was sufficient under the current method. The inclusion 
was discretized into 27 smaller cubes {a = 27) within which the eigen- 
stresses and eigenvelocities were approximated as constants. A total 
of 23^ — 1 fourier terms were used in the expansion of the periodic 
parts of the field variables. 

It should be noted that all the components of the 4 effective tensors 
are real. This is due to the fact that the unit cell is symmetric in all 
three orthogonal directions. The essential relation between the cou- 
pling tensors ((S^)^ = S^) is verified from the present calculation. The 
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Table 1: Effective properties calculated for Q = 1.5 on branch 



effective compliance tensor D has both major and minor symmetries 
{Dijki = Dijik = Djiki = Dkiij) due to the symmetry of the problem. 
In general we should have Dij^i = Dl^^-y Due to the symmetry of the 
unit cell, the effective density matrix is diagonal and indicates that 
there is no coupling between the average momentum in a particular 
direction and the average velocity in an orthogonal direction. The 
components of all the four effective tensors indicate that there is a 
distinction between the x direction and the other two directions by 
virtue of the fact that the propagation direction is x. On the other 
hand the other two directions, y and z, are indistinguishable for a 
wave propagating in the x direction due to the symmetry of the unit 
cell. This is seen to be true from the calculated components of the 
effective tensors (eg. ^2222 = ^3333 ¥= -Dim etc.) 

4 Conclusions 

A method for the homogenization of 3-D periodic clastic composites 
is presented. It allows for the evaluation of the overall dynamic mate- 
rial constitutive tensors relating the averaged dynamic field variable 
tensors of velocity, strain, stress, and linear momentum. The method 
doesn't require the apriori calculation of the solution to the elasto- 
dynamic problem. Rather, the dispersion relation and the pointwise 
solution of the composite emerge from the current method. The cou- 
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pled form of the constitutive relation as proposed by Milton and Willis 
emerge natTirally from the present method. We have also shown that 
the matrices corresponding to the effective comliance and density ten- 
sors are hermitian and that the coupling tensors are hermitian transose 
of each other. Finally we have presented corroboration of the theoret- 
ical proofs in the form of a solved example of a 3-D periodic composite 
composed of elastic cubes periodically distributed in an elastic matrix. 
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A Appendix 

For isotropic case, 

Cfjki = >^'^SijSki + A^ikSji + SuSjk] (57) 
The left hand side of Eq. (fT3|) can be expanded as, 



-C • C° : (C ® u)sym + wV°u = Hijiij = (58) 

[-(AO + ;,o)c.o + (^y - A')^.,] % 



To invert the tensor H in Eq. (|58j) we seek a tensor J such that, 

Jmi • Hjjr' — -^mjr — ^mj (59) 

Expressing J in terms of C (8) C ^-iid l'^ we have, 

J = a(C®C) + 61^ (60) 
Contracting with H we get the following equation, 

[-a(AO + /)C2 + a(^V° - M°C') - + /)] C™0+K^V°- A')^m, 
The constants, therefore, are. 



[a;2p° - M°C^] - (AO + 2/xO)C2] 

6 



(61) 



Now Eq. (|13p can be written as. 



u(0 = *(0-U(0 + ©(0:5](^) (62) 

where 

= [aCiO + ft<5*j] (63) 

and 
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Since gets contracted with the symmetric tensor XI, it can be 
symmetrized about its last two indices. We make the following trans- 
formation, 



(65) 



[^ijCp + ^ipQ] 



2 [AO + ^0 



0^2] CiCiCp 



0>2 ["^uCp + ^ipCj] 



The left hand side of Eq. (|14p can be expanded as, 



-C : [C «) (C • (T)]sym + UJ p a- = Fijkiaki 
.1 



(< 



2^0s 



where we define the following tensors. 



dijli = [Ci^jkCi + Ci^jiCk + Cj^ikCi + Cj^iiCk 



(67) 



9ijkl ~ ^ijCkCl'i 9ijll — CiCj^kl', Qijli — CiCjCkCl 

and note the following relations. 
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mnij^ijkl 
(1) J2) 



(1) 



(1) 



(2) 

(2) (3) 

(3) d) 



Hmniji'ijkl 
(2) 

5. 



.(3) 



(3) 



mmjdi^kl 

(4) d) 

mniji'ijkl 
Hmnij^ijkl 



ymnfci 
^ ^mnkl 
C^^mn^kl 

"^ymnkl 

'^ilmnkl 
^2 (3) 
^ ^mnkl 

^ ^mnkl 



(69) 



mnkl 



To invert the tensor F in Eq. (j66p we seek a tensor G such that 



G : F 



I 4s 



Expressing G in terms of g^*^ and l^** we have, 

G = ag(i) + 6g(2) + cg(3) + dg(4) + el^' 
Contracting with F we get the following equation, 



(70) 



(71) 



[a(c.V° - /C) - e/] g(^) + [6(a;V° - (A° + 2/)C2) - eA"] g 
+ [ca;V°] g^'^ + [-a(2A0 + 2/) - c(3A0 + 2/) - d{io^ - (A^ + 2/)C2)] g 

Solving the above equation we get the following values for the 
constants 
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c = 
d= - 



M 

>^ 

2^0 [AO + ^0] 



Now Eq. ()14p can be written as, 

a(0 = *(0-U(0+r(^):S(0 

where, 



(72) 



(73) 



il/ . . 



+ 
and 



2/i0(A0 + /) 



[^^2p° - /uOC^][wV° - (A° + 2pO)C^] 



CiCjCp 



(74) 
(75) 



AO 



a;2p0 _ (AO + 2/^0)^2 



a;2p0 _ ^0^2 



r 



ijkl 



+ 



+ 



/iO 1 

~^n\ ?772oi''''^ifc^' + Ci^jlCk + Of^ifcO + CjklCk} 



(76) 
(77) 



AO 



tj2p0 _ (AO + 2^0)C^ 

2/x0[A0 + /i0] 



^ijCkCl 



[w2p0_^0^2][^2^0_(A0 + 2/z0)CV'^" ' 2 



CiCjCkCl + T; i^ik^jl + ^il^jk] 



The components of tensors 0, T, and ^ can be expressed in 
terms of the longitudinal wave speed ci = (A^ + 2pfi)/ and the 
shear wave speed ci = ^y/jP/fP, 



r2 
^2 



[^2 _ c2C2] [a;2 - cfC 



,2>21 ''^'=1 



[^2 _ e2c2 



(78) 
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(79) 



ijkl 



4 1 



+ 



{Ci^jkCl + Ci^jlCk + Cj^ikCl + Cj^ilCk} 



(80) 



a;2 - dC 



,2/-2 



^ijCkCl + 



[a;2-ciC2][a;2-c2c2]— ' 2 



IJP 



-up 



24(4 



+ 



cf - 2c: 



|o;2_c2c2][a;2-cK 

2 



2 /-Si 



CiCjCp 



(81) 



To facilitate solution we convert the tensorial equations into their 
matrix form in the main text which are then used to prove some es- 
sential behaviour of the tensors. These proofs rely on transforming 
the tensors into forms where their inter-relationships are more trans- 
parent. To that end we transform the tensor ^ to, 



where D° is the compliance tensor given by 



(82) 



d: 







-A' 







1 



'""^^ 2/x0(3A0 + 2/xU)"""''^ ■ Afj,' 
The final expression for 4f is given below, 



mnp 



nO ..xCr. 
~2 



(84) 



2 [cl - 4] 



- 4e] [c.^ - 4e] 



CmCnCp 



+ 



- 4e 

If we now define as. 



mp\ 



®mnp P ®mnp 



(85) 
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we find from Eqs. (|65l [Hi]) that 



We also define ^ as, 



and f as, 



* mnp ^pmn 



(86) 



(87) 



+ 



4(a;2-ciC' 
(cf - 2ci 



■{CmSnkCl + CmSnlCk + Cn^mkQ + Cn^mlCk} 



2ci(3c?-4ci 
1 



We notice that transformed as above, T has the symmetries, 

^mnkl ^nmkl ^mnlk ^klmn 
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